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SECTION  I 
INTRODUCTION 


Electromagnetic  radiation  in  the  presence  of  stratified  media  is  a 
fundamental  problem  in  geophysical  sciences,  printed  circuits,  microstrip 
devices,  etc.  [1,2].  It  has  been  a  subject  of  extensive  research  for 
eighty  years  beginning  with  the  well-known  Sommerfeld  integral.  However, 
this  is  also  an  area  noted  for  its  lack  of  nimerical  data  and  practical, 
generalized  solutions.  The  electromagnetic  problem  is  complicated  because 
of  a  variety  of  physical  parameters  such  as  the  source  location,  the  layer 
thicknesses,  and  dielectric  constants  in  individual  layers.  The  solutions 
are  often  complex  and  are,  as  a  rule,  only  applicable  to  particular  cases. 

With  the  advent  of  high-speed  digital  computers,  it  is  increasingly 
more  desirable  to  seek  numerical  solutions  that  are  simple  to  use  and  are 
applicable  to  a  broad  spectrum  of  problems.  The  generalized  formulation  in 
terms  of  reflection  coefficients  and  one-dimensional  integrals  in  k 

P 

appears  to  be  a  very  suitable  representation  for  numerical  computations 
[3,4].  The  primary  difficulty  in  this  approach  is  in  the  computation  of 
the  integrals  which  involve  Hankel  functions. 

A  highly  efficient  technique  to  compute  these  integrals  is  by  way  of 
Fast  Fourier  Transform  (FFT)  [  5-8  I.  The  integrand  is  organized  into  a 
product  of  two  functions,  one  of  which  includes  a  Bessel  function.  The 
function  without  the  Bessel  function  is  expanded  into  a  discrete  Fourier 
series  via  FFT.  The  exponential  terms  of  each  of  the  resulting  discrete 
Fourier  series  containing  kp  are  then  combined  with  the  remaining  integrand 
for  closed-form  integration.  There  are  three  difficulties  involved  in  this 
approach.  Firstly,  the  integral  after  FFT  expansion  must  be  suitable  for 
evaluation  in  closed  form.  Secondly,  the  resulting  closed-form  integral 
must  be  easy  to  compute  numerically  since  there  may  be  thousands  of  them  in 
the  Fourier  series  expansion.  Thirdly,  proper  decay  of  the  integrand  is 
required  for  the  FFT  expansion. 

In  this  report  this  approach  is  expanded  to  cover  the  general  problem 
of  radiation  and  scattering  in  the  presence  of  stratified  media.  The 
integrand  is,  in  most  practical  cases,  a  product  of  two  functions;  an  even 
or  odd  function  of  k  and  a  Hankel  function  (or  its  derivative)  of  the 


first  kind.  It  is  demonstrated  in  this  report  that  a  general  algorithm  to 
compute  this  class  of  problem  is  feasible. 
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SECTION  II 

FORMULATION  OF  THE  PROBLEM 


The  fields  in  the  &-th  layer  due  to  a  Hertzian  dipole  in  a  stratified 
aedim  as  shown  in  Figure  1  can  be  expressed  in  the  following  general  form 
[4  ] 


dk  [AAe^V  +  Bfeejk^2)  Hn  V)Cn4) 


^  [  -A*e-jk**Z  +  B£ejk^2]  H^1)(kpp)Cn(0) 


dk  f  C.e  Jk*z  +  Doejk^z  1  H(1)(k 

p“k2T  1  n  P  n 

P 

(2) 


dk  [  -A*e  Jk^2  +  B^fcz2  ]  H(1^k  p)C’  (d>) 

p  .  2  1  *  *■  Jnpn 

k  P 
P 


dk  [  C«e"jkAz2  +  Dzejk*z2  ]  H(1),(k  p)S  (<p) 

n  c  1  *  K  ‘  n  p  n 


»«*  -J 


dk  [  C2e"JkWz  +  DLeiklz*  ]  H(1)(k  p)S  (*) 
p  *  *>  npn 
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H*p  « 
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dk  -7——  [  -C^e"^**2  +  D  ejk£zZ]  H^1)(k^p)S^(<^> 


p  K 


n  p  n 


dkp  “kV^  1  A^e"jk£zZ  +  BjiejkizZ]  H^1)(kpp)c;(^>) 
P 

(5) 


dk 


P  ~j~ 

kpP 


[  -Cie“jk£zZ  +  Diejk*zZ  ]  H^1)  o)S^(*) 


dk 


[  Aje-^**2  +  Biejk*zZ)  H(1)'(k  p)C  (<t>) 


n  p  n 

(6) 


where  denotes  the  n-th  order  Hankel  function  of  the  first  kind  and 

n  n 

its  derivative.  Sn,  Cn  and  n  are  determined  by  the  dipole  excitation 

involved.  S'  and  C'  are  derivatives  of  S  and  C  .  The  e  JU)t  term  is 
n  n  n  n 

suppressed  and  the  cylindrical  coordinates  are  adopted ,  in  which 


k 


2 

i 


k 


(7) 


for  the  i^-th  layer.  In  Equation  (7),  k  is  independent  of  the  layer  and 

kl  ”  “  *  P 

The  four  constants,  A^,  B^,  Cjg,,  and  D jj,  are  determined  by  enforcing  the 
boundary  conditions  at  each  layer  interface.  For  example,  the  continuity 
of  displacement  flux  5  normal  to  the  interface  at  z  ■  d £  requires 


et  (AAe“Jk*zd£  +  Bjtejk4zdA  ) 

-  e£_1(AjUle“jk(i-l)zdi  +  BA_1eJk(A-i)zdA) 


(8) 
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There  are  a  total  of  four  equations  similar  to  Equation  (8)  at  each 
interface  which  constitute  the  necessary  and  sufficient  boundary 
conditions . 

The  quantities  n,  Sn,  Cn>  Aq,  Boi  C0,  D0  are  determined  by  the  dipole 
excitation.  First,  the  fields  due  to  a  Hertzian  dipole  in  the  absence  of 
stratification  are  expressed  in  integral  forms  similar  to  Equations  (1-6). 
These  constants  and  coefficients  are  then  determined  by  a  comparison 
between  Equations  (l)-(6)  for  A  *  0  and  the  fields  for  this  particular 
dipole  source  in  the  absence  of  stratification.  The  details  are  included 
in  References  3  and  4,  and  are  included  in  Appendix  A  for  convenience. 
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SECTION  III 

FFT  METHOD  AS  APPLIED  TO  A  SPECIAL  CASE 


Before  presenting  the  application  of  FFT  method  to  the  general 
integrals,  the  FFT  method  that  was  used  [5],  [6]  to  evaluate  a  specific 
integral  will  be  discussed.  The  integral  is  given  as, 


•  L 


E(k  )  (k  p)  dk 
pip  p 


-/ 


E(k  )  J,(k  p)  dk 
P  1  P  p 


where  E  denotes  an  even  function.  We  first  define 


(8) 


g(kp)  -  2E(kp)  exp(vRkp)  . 


(9) 


The  Discrete  Fourier  Transform  (DFT)  of  g  can  then  be  expressed  as, 

G(f )  -  J  g(kj  exp  (-2*jfkJ  dk„  ,  (10) 


g(kp)  exp  (-2itjfkp)  dkp  , 


where 


f  -  n/(Ndk  ). 
P 


We  can  now  approximate  g(kp)  by  a  DFT  expansion  as  follows: 
,  (  N/2)  -1 

g(kp)%N kr  ^  G(f)exp(j2Trfkp) 

n*-N/2 


(11) 


for  0  <  k  <[(N/2)-  1  ]A  kp 

By  substituting  Equation  (11)  into  Equation  (8)  and  recognizing  that 
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/: 


e“VkP  J. (k  p)dk  -  —  [1  - 
1  P  p  p 


(v*  +  p>)1/2 


]  . 


(12) 


where  v  »  vR  +  jVj,  Re(v  +  jp)<  0  and  we  have. 


IE  1  ~  1  r(  n  X  _L 

Hi  NAkp  /  ^NAkp^  p 

n*-N/2 

U  "  (V  NAk?  )/l(VR  "  J2^/(NAkp))2  +  p2)1/2}  . 


(13) 


Thus  Che  integral  in  Equation  (8)  is  transformed  into  a  finite  series 
expansion  in  Equation  (13),  and  can  be  readily  computed  on  a  digital 
computer . 

Several  fine  points  in  this  method  need  careful  and  proper  handling. 

An  arbitrarily  chosen  nimber,  vR  is  required  in  Equation  (12)  to  be  larger 
than  zero.  However,  VR  must  also  be  smaller  than  z  so  that  g(k^)  as  defined 
in  Equation  (9)  vanishes  as  kp  approaches  infinity.  In  the  problem  treated 
in  Reference  5,  a  term  exp  (jk^z  ♦  vR  kp)  is  contained  in  g(kp).  In  order 
to  maintain  a  finite  g(kp),  vR  ■  Z/2  is  chosen  so  that  g(k  )  vanishes  as  k 
approaches  infinity.  Another  concern  is  that  a  real  would  lead  to 
singularities  in  the-integrand .  This  difficulty  can  be  overcome  by  adding 
a  small  dielectric  loss  to  remove  the  singularity  from  the  real  axis. 

A  key  feature  of  this  approach  is  to  transform  the  interval  of 

integration  (—*,«)  into  (o,  “)  so  that  g(k  )  nearly  vanishes  outside  (o, 

N  P 

(j  -1)  Akp)  for  a  properly  chosen  VR.  After  this  is  accomplished,  closed- 
form  integration  similar  to  that  of  Equation  (12)  is  needed  to  handle  the 
remaining  integral  containing  a  Bessel  function.  In  essence,  the  FFT 
method  transforms  an  integral  containing  a  Bessel  function  to  an 
approximate  summation  which  is  void  of  special  functions  and  therefore  can 
be  readily  computed.  In  the  next  section,  the  FFT  method  will  be  shown  to 
be  applicable  to  all  the  integrals  in  Equations  (1-6).  This  can  be 
established  by  showing  that  for  these  integrals  the  interval  of  integration 
can  be  changed  from  (  ~°°,  “)  to  (o,  »  )  and  simple,  closed-form  integration 
can  be  carried  out  for  the  remaining  integral  containing  a  Bessel  function. 
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SECTION  IV 

APPLICATION  OF  FFT  METHOD  TO  THE  GENERAL  INTEGRALS 


For  dipole  excitation,  the  order  of  the  Hankel  functon,  n,  is  either  0 
or  1.  The  coefficients  A^,  B^,  C  ,  and  are  all  formulated  in  matrices 
and  can  be  readily  evaluated  on  a  computer.  The  evaluation  of  the 
integrals  involving  Hankel  functions  can  be  handled  systematically  by  the 
FFT  method  just  described. 

A  key  to  the  success  of  the  FFT  method  is  the  term  exp(jk^  z  + 
vr  kp) ,  which  causes  g(kp)  to  vanish  for  a  properly  chosen  positive  vr,  as 
k^  approaches  infinity.  As  can  be  seen  in  the  Appendix,  terms  of  this  type 
exist  in  the  general  form  of  the  integral.  That  a  suitable  vr  exists  and 
the  integrand  would  vanish  as  k^  approaches  inf inity  can  be  argued  on  the 
basis  of  the  existence  of  the  solutions.  It  can  also  be  proved,  separately 
for  the  cases  of  a  dipole  below  and  above  the  field  point,  that  a  suitable 
vR  exists . 

The  next  hindrance  to  the  FFT  method  is  in  the  closed-form  evaluation 
of  the  remaining  integral  containing  the  Bessel  function.  It  is  shown  in 
the  following  that  a  closed-form  integration  similar  to  Equation  (12)  does 
exist  for  each  integral  that  may  be  encountered. 

The  integrals  involved  in  Equations  (1-6)  can  be  divided  into  the 
following  eight  types,  each  of  which  can  be  handled  in  a  manner  similar  to 
that  of  Equation  (8)  discussed  in  the  preceding  section.  E  and  0  denote 
"even"  and  "odd"  functions  in  the  following  discussion. 


(A)  I* 

Ho 

For  the  integral 


E(k  )H(1)(k  p)dk 

p  o  p  p 


E(k  )H  (k  p)dk  , 

pop  p 


(14) 
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where  N0  is  a  Bessel  function  of  the  second  kind  of  order  0,  the  following 
closed-fora  integral  can  be  eaployed  , 


e~vkp  N  (k  p)dk  -  - — 


24.„2vi/2 


V+(vZ+pQ 


P  P  r(v2+p2)1/2 


where 


Re  v  >  I  p  *  0  . 

ID 


(b)  i: 


For  the  integral 


v  /. 


0(k  )H(1)(k  p)dk 

P  o  p  p 


«2  /  0  (k  )  J  (k  p )  dk  , 

/  o  POP  P 


the  following  closed-form  integral  can  be  employed 


e  Vkp  J  (k  p)dk  -  1  -  - - - r-py~  , 

°  °  P  (v2+p2)1/2 


where 


Re  [v  +  jp]  >  0 


(C)  xh' 
"o 


For  the  integral 


'  /. 


E(k  )H(1)'(k  p)dk 

P  o  p  p 


E(kp)H'1)(kp»)dkp 
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/; 


EOcUrkpJdk  , 

Pip  p 


(18) 


the  closed-form  integral  of  Equation  (12)  can  be  used. 


(D)  I® 
Ho 


Since 


I0  = 


/: 

7- 


0(k  )H(1)' (k  p)dk 
p  o  P  p 


0(k  )H(1)(k  p)dk 

pi  p  p 


2 '  V 

the  discussion  of  1°  later  in  (F)  will  be  applicable. 


(19) 


Hi 


(E)  I 


«1 


This  type  of  integral  has  been  discussed  in  Section  III. 


(F)  I 


Hi 


For  the  integral 


'  /- 


0(k  )  H(1)(k  P)dk 

p  1  p  p 


/ 


-  2j  |  0(k  )N. (k  p)dk 

pip  p 


.2j  [^o(kp)No(kpp)]  fc 


0  J  P 

J  o 


(k  P)0' (k  )dk  , 
op  P  P 

(20) 
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where  0'(kp)  is  the  derivative  of  0(kp),  the  following  integral  can  be 
employed 


w)dk» 


_ zZ - tn  inSyHe }£!1 

.<v2+(>2>l/2 


where 


Rev  >  Ithp  *  0  . 


(21) 


In  the  last  step  of  Equation  (20),  the  first  term  can  be  readily 
evaluated  by  substitution.  However,  this  term  generally  vanishes  because 
of  the  unique  nature  of  N0  at  o  and  00  . 


(G)  I*. 


Since 


*  ■  /. 

/00 

-CO 

7- 


E(k  )H(1)'(k  p)dk 
Pi  P  P 


E(k  )H(1)(k  p)dk 
p  o  p  p 


r^r  E(k  )H(1)(k  p)dk 
kpp  pi  p  p 


E  0 

previous  discussions  on  I„  and  Iu  are  applicable. 

Ho  M1 


0(k  )H.(1)’(k  p)dk 
Pi  P  P 

0(k  )H<X)(k  p)dk 

p  o  p  p 

7~  0(k  )H1(1)(k  p)dk  , 
k  p  pi  p  p 


(22) 


(23) 
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SECTION  V 

NUMERICAL  ANALYSIS 


The  preceding  discussions  demonstrated  that  a  general  numerical 
approach  for  the  problem  of  dipole  radiation  in  stratified  media  is 
feasible.  In  fact,  a  computer  algorithm  can  be  developed  for  the  general 
case  of  arbitrary  dipole  configurations,  geometries  of  stratification,  and 
dielectric  and  ferromagnetic  properties.  The  well-organized  integrals  in 
Equations  (1-6)  can  be  readily  computed  by  the  FFT  method  after  the 
parameters  are  specified.  The  matrix  formulation  for  the  coefficients  A^, 
B^,  C^,  and  can  be  easily  computed  for  any  number  of  layers,  dielectric 
constants,  layer  thicknesses,  etc.  Note  also  that  these  four  coefficients 
are  all  even  functions  of  kp  .  Whether  the  integrands  devoid  of  the  Bessel 
functions  in  Equations  (1-6)  are  odd  or  even  can  also  be  readily 
determined . 

Several  delicate  points  in  this  method  should  be  carefully  observed. 
The  choice  of  N,  Ak^  and  must  be  properly  made  to  insure  that  the  DFT 

expansion  is  valid.  As  will  be  shown  in  the  next  section,  good  convergence 
and  accuracy  were  achieved  in  all  the  examples  tested.  Obviously,  the  fact 
that  the  basic  FFT  techniques  are  well  established  has  a  great  deal  to  do 
with  the  success  of  the  present  method. 

A  major  limitation  of  this  method  is  its  deteriorated  convergence  in 
computing  fields  at  locations  with  large  P  .  In  order  that  g(kp  )  be 
sufficiently  near  zero  outside  (-NAkp/2,  NAk^/2)  for  the  FFT  expansion,  Akp 
must  be  chosen  to  be  much  smaller  than  1/p  ;  and  N  must  be  chosen 
sufficiently  large.  The  limited  central  memory  and  speed  of  modern  digital 
computers  restrict  the  FFT  expansion  to  finite  distances  of  p  .  Depending 
on  the  nature  of  the  function  g(k^  ),  accuracy  of  the  computation 
deteriorates  when  p  is  larger  than  several  to  several  hundred  wavelengths. 
Extension  of  this  method  to  cases  of  large  p  is  possible  and,  of  course, 
highly  desirable. 

In  this  research,  no  attempt  was  made  to  develop  a  generalized 
computer  program.  Instead,  basic  techniques  and  subroutines  were  developed 
and  tested  to  demonstrate  that  a  generalized  algorithm  is  indeed  feasible. 
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SECTION  VI 
NUMERICAL  TESTING 

Several  methods  were  employed  in  testing  the  accuracy,  efficiency  and 
versatility  of  the  FFT  method  for  problems  involving  stratified  media.  The 
results  were  checked  against  known  data  in  the  literature,  against  data 
obtained  for  limiting  cases,  and  by  numerical  convergence  with  variations 

in  N  and  Ak  .  All  the  three  field  components,  Ez,  E  ,  and  E  ,  were 

p  p  _  p 

computed  for  the  more  involved  cases  of  horizontal  electric  and  magnetic 

Hertzian  dipoles.  The  excellent  numerical  results  demonstrated  that  the 
present  method  is  indeed  a  powerful  and  general  computing  tool  for 
electromagnetic  problems  in  stratified  media. 

The  first  case  tested  was  for  Ez  due  to  a  horizontal  electric  dipole 
in  media  of  two  and  three  layers,  which  had  previously  been  computed  by  the 
FFT  method  [5,6] .  Excellent  agreement  between  the  present  analysis  and  FFT 
results  was  observed. 

Ez  due  to  a  horizontal  magnetic  dipole  is  similar  to  that  due  to  a 
horizontal  electric  dipole  because  they  both  have  the  same  type  of 

j? 

integral,  Iu  of  Equation  (8).  Figure  2  shows  a  case  that  has  been 
"1 

extensively  studied  for  t  *  0  and  £2  *  80eo  by  Kong  and  his  colleagues  [9] 

-[11].  Their  earlier  results  in  References  9  and  10  were  inadequate  and 

were  later  discussed  [11]  in  conjunction  with  the  newly  obtained  correct 

data.  It  must  be  emphasized  that  the  integral  involved  is  of  the  same  type 

as  that  computed  by  FFT  method  in  their  earlier  work.  Even  though  some 

convergence  problem  may  exist  since  the  source  and  field  points  have  the 

same  z  coordinates,  their  earlier  FFT  algorithm  should  have  readily  yielded 

limiting-case  data  if  the  general  and  universal  nature  of  this  method  had 

been  recognized  and  taken  advantage  of.  This  is  illustrated  in  Figure  3, 

in  which  a  comparison  is  made  for  Ez  with  various  values  of  t  to  compare 

with  the  data  for  t  *  0  in  the  literature.  The  present  computation  was 

actually  made  for  the  case  in  which  the  layer  with  water,  for  which  eo~80e  , 

i  o 

was  approximated  by  a  perfect  conductor  as  often  done  in  experiments. 

As  can  be  seen,  the  computed  results  converge  nicely  to  the  t  ■  0 

case  when  t  is  reduced  from  0.1  A  to  0.02  A  .  The  ripples  in  the  curve 

0  0 

corresponding  to  t  *  0.02 A  are  indications  of  poor  numerical  convergence 
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Figure  2.  A  horizontal  magnetic  Hertzian  dipole  in  a 
three-layer  medium. 


800 


t  ^  0  approach  those  for 


which  will  be  discussed  later  in  greater  detail.  The  deterioration  of 
convergence  for  small  t  is  due  to  the  reduced  distance  in  z  coordinate 
between  the  source  dipole  and  the  field  point,  as  was  explained  in  an 
earlier  section. 

Because  of  a  lack  of  existing  data  that  could  be  used  to  check  against 
the  present  computation  involving  horizontal  magnetic  dipoles,  we  resorted 
to  cases  for  which  independent  image-theory  calculations  can  be  made.  For 
example,  a  horizontal  magnetic  Hertzian  dipole  in  a  homogeneous  medium  over 
a  ground  plane  can  be  analyzed  easily  by  image  theory.  These  data  can  be 
used  to  compare  with  those  computed  by  the  FFT  method.  Figure  4  shows  a 
comparison  in  the  amplitude  and  phase  of  Ep  for  a  horizontal  magnetic 
dipole  in  a  medium  with  e  ■  2  e0  over  a  conducting  plane.  Figure  5  shows 
a  comparison  for  E.  in  this  case.  An  accuracy  of  0.1  percent  or  better  is 
achieved  in  these  cases . 

The  convergence  test  is  another  useful  check  in  the  present  FFT 

method.  Figure  6  shows  how  convergence  is  reached  by  increasing  N  in  the 

computation.  It  was  observed  that  numerical  convergence  is  consistently 

accompanied  by  the  smoothing  of  the  ripples  in  the  E  vs.  p  plots. 

Convergence  can  also  be  checked  by  varying  the  value  of  Ak  chosen  in  the 

P 

FFT  computation. 

All  types  of  integrals  involved  in  Equations  (1-6)  and  discussed  in 
Section  IV  were  numerically  tested  successfully.  Thus,  the  applicability 
of  the  FFT  method  to  the  general  electromagnetic  radiation  in  stratified 
media  has  been  clearly  demonstrated. 
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Figure  4.  Comparison  between  the  exact  image-theory  calculation  and  the 

present  FFT  method  for  E  due  to  a  horizontal  magnetic  Hertzian 
dipole  in  a  medium  of  2e  bounded  by  an  infinite  conducting  plane 
one  free-space  wavelengtR  below. 
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SECTION  VII 
CONCLUDING  REMARKS 


It  is  established  in  this  report  that  a  general  and  universal  method 
exists  for  numerical  analyses  of  the  general  electromagnetic  problem  in 
stratified  media.  The  problem  can  be  formulated  such  that  a  computer 
algorithm  can  be  written  for  the  general  problem  of  arbitrary  excitation, 
layer  thickness,  dielectric  constants,  etc.  The  integration  is  by  the  FFT 
method  which  transforms  the  integrals  involving  Bessel  functions  into  a 
summation  of  simpler  functions.  Excellent  numerical  results  were  obtained 
for  many  test  cases  some  of  which  are  presented  in  this  report. 

There  are,  however,  two  shortcomings  in  the  present  FFT  method. 

Because  of  the  limited  central  memory  and  execution  speed,  computation  is 
limited  to  field  points  having  a  radial  (p)  coordinate  of  several  to 
several  hundred  wavelengths,  depending  on  the  nature  of  the  particular 
integrand  involved.  In  addition,  convergence  deteriorates  as  the  distance 
(in  z-coordinate )  between  the  source  and  field  points  decreases. 

For  future  research,  it  is  recoianended  that  a  user-oriented  general 
computer  algorithm  be  developed  on  the  basis  of  the  findings  in  this 
report.  Furthermore,  the  two  shortcomings  in  this  FFT  method  should  be 
examined  for  improvements. 
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APPENDIX  A 


CONSTANTS  AND  COEFFICIENTS  FOR  THE  kp- space  INTEGRAL 
EXPRESSIONS  OF  FIELDS  IN  STRATIFIED  MEDIA 


Equations  (1-6)  are  valid  for  the  general  problem  of  dipole  radiation 
in  stratified  media.  Although  the  application  of  these  results  to  a 
specific  case  is  straightforward,  it  is  convenient  to  sumnarize  all  the 
relations  needed  in  the  determination  of  all  the  constants  and 
coefficients,  which  are  described  in  detail  in  Reference  4. 

At  each  medium  interface  the  following  boundary  conditions  exist 

e^A^e  JkJi.zdi  +  B£e-*kJtzdS.) 
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(A-4) 


An  arbitrary  Hertzian  dipole  can  be  considered  to  be  the  superposition 
of  four  basic  Hertzian  dipoles,  vertical  electric,  horizontal  electric, 
vertical  magnetic,  and  horizontal  magnetic.  The  quantities  n,  Cn,  Sn,  Aq, 
Bo*  CQ,  ®o»  are  determined  by  comparing  Equations  (1-6)  with  the  fields  of 
the  specific  dipole  under  consideration  with  the  fields  due  to  this 
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particular  dipole  in  the  absence  of  stratification  (infinite  region).  The 
latter  is  often  referred  to  as  the  "primary"  contribution  in  the  general 
strati fied-media  solution  and  is  summarized  as  follows 


(1)  Vertical  electric  dipole  (VED) 
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and  IA  is  the  electric  dipole  moment.  Note  here  that  0  denotes  layer 
index,  not  "free  space",  in  kQ,  Pq* 


(2)  Horizontal  electric  dipole  (HED) 
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(3)  Vertical  magnetic  dipole  (VMD) 
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and  1A  is  the  magnetic  dipole  moment. 


(4)  Horizontal  magnetic  dipole  (HMD) 
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Because  of  the  functional  difference  between  z  >  0  and  z 
integrands  of  the  fields  described  above,  the  coefficients  Aq, 
Do  are  denoted  by  Aq+,  Bq+,  Cq+,  and  Dq+  for  z  L  0  and  Aq~,  Bq 
Dq-  for  z  <_  0. 
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The  fields  in  region  0  for  the  four  types  of  dipoles  cen  be  determined 
by  the  following  relations 


(1)  VED 
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for  z  >  0  (A- 19) 


for  z  <  0  (A-20) 


where  Avfid  and  Bved  characterize  contributions  due  to  the  stratified  medium 
and  are  to  be  determined  by  the  boundary  conditions. 


(2)  HED 
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for  z  1  0  (A-22) 


where  A^ed,  ®he<j»  Ched*  and  Dhed  8re  to  Ee  determined  hy  the  boundary 
conditions . 
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(3)  VMD 
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where  and  are  to  be  determined  by  the  boundary  conditions. 

(4)  HMD 
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for  z  <  0 


(A-26) 


where  ^d*  Bhmd’  Chmd  ’  and  Dhmd  are  to  be  deter“ined  by  the  boundary 
conditions . 

It  is  convenient  to  reorganize  Equations  (A-l)-(A-4)  in  either 
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or 
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and 
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The  wave  amplitudes  in  region  0  for  the  four  types  of 
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(4)  for  horizontal  magnetic  dipole 
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